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ABSTRACT 



Aims. The main goal of this analysis is to present a new method to estimate the physical properties of diffuse cloud of 
atomic hydrogen observed at high Galactic latitude. 

Methods. This method, based on a comparison of the observations with fractional Brownian motion simulations, uses 
the statistical properties of the integrated emission, centroid velocity and line width to constrain the physical properties 
of the 3D density and velocity fields, as well as the average temperature of Hi. 

Results. We applied this method to interpret 21 cm observations obtained with the Green Bank Telescope of a very 
diffuse HI cloud at high Galactic latitude located in Firback North 1. We first show that the observations cannot be 
reproduced solely by highly-turbulent CNM type gas and that there is a significant contribution of thermal broadening 
to the line width observed. To reproduce the profiles one needs to invoke two components with different average 
temperature and filling factor. We established that, in this very diffuse part of the ISM, 2/3 of the column density is 
made of WNM and 1/3 of thermally unstable gas (T ~ 2600 K). The WNM gas is mildly supersonic (< M >~ 1) and 
the unstable phase is definitely sub-sonic (< M >~ 0.3). The density contrast (i.e., the standard deviation relative 
to the mean of density distribution) of both components is close to 0.8. The filling factor of the WNM is 10 times 
higher that of the unstable gas, which has a density structure closer to what would be expected for CNM gas. This 



field contains a signature of CNM type gas at a very low level (Nh 
convergent flow of WNM gas. 



3 x 10 ) which could have been formed by a 
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1. Introduction 

Neutral atomic hydrogen (Hi) is the most abundant phase 
of the interstellar medium (ISM) and it plays a key role 
in the cycle of interstellar matter from hot and diffuse 
plasma to the formation of stars. The 21 cm transition 
has been used sinc e the 1950s (jEwen fc Purcelll 119511 
iMuller fc OortJ . il 95ll ) to study the properties of Galactic Hi 
but some of its basic properties are still unknow n. The gen- 
eral p icture of i nterstellar Hl first modele d bv lField et al.l 
(fl969l) (see also iWolfire eTail (U995L \2tM) ). is of spatially 
confined cold structures (CNM - Cold Neutral Medium, 
T sa 100 K) surrounded by a more diffuse and warmer 
phase (WNM - Warm Neutral Medium, T » 8000 K). The 
heating and cooling mechanisms in the diffuse ISM are such 
that these two phases can co-exist in pressure equilibrium. 
On the other hand the exact physical conditions in which 
the thermal transition from WNM to CNM occurs - most 
probably related to a local increase of the density and to 
the temperature dependence of the cooling function - are 
still unclear. Turbulence plays a dominant role in the kine- 
matics and density structure of Hi: both phases (CNM and 
WNM) have a self-similar density and velocity structure 
and an energy spectrum compa tible with a turbulent flow 
(|Miville-Deschenes et all . l2003al) . 
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With the advent of more sensitive and higher angu- 
lar resolution observations, new challenges have come up 
(e.g., tiny scale atomic structures - TSAS -, self- absorption 
structures and thermally unstable gas). An important ad- 
vance came re cently from the M i llenniu m Arecibo Survey 
conducted by iHeiles fc Trolandl (|2003al fbh. These authors 
claimed that a significant fraction (~ 30%) of the Hi is 
in a thermally unstable regime. This observational evi- 
dence seems to be in accord with recent numerical simu- 
lations devoted to the stud y of the Hi thermal instability 
(| Audit fe Hennebelld . l2005h . 

Here we contribute to this effort by analyzing 21 cm 
observations of a very diffuse high-latitude region of the 
Galactic ISM. The main advantage of studying the proper- 
ties of weak Hi feature is that no modeling of radiative 
transfer is required to interpret the brightness observed 
which can be converted directly in Hi column density. 
Furthermore the 21 cm emission of such high latitude re- 
gions do not suffer from velocity crowding which facilitates 
the decomposition of the spectra. Our goal is to study the 
dynamical and thermo-dynamical properties of the neutral 
gas in a diffuse region of the solar neighborhood. 

Unlike previous analysis of 21 cm spectra we do not use 
a combination of emission and absorption observations to 
study the thermal and kinematical properties of the Hi but 
instead we use the statistical properties of the emission on 
the plane of the sky to infer the three-dimensional (3D) 
properties of the gas. In particular we exploit the idea that 
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variations of the value of the velocity centroid on the plane 
of the sky can be used to constrain the amplitude of the 
turbulent motions on the line of sight and therefore put 
limits on the average thermal broadening over the field. 

In §[2] we describe the 21 cm observations used for this 
analysis and some quantities computed from them that will 
be used further to infer the Hi physical properties. Then we 
describe (§[3]) the method used to simulate realistic 21 cm 
observations using fractional Brownian motion and to test 
the hypothesis of an isothermal and turbulent Hi. In §0] we 
estimate the Hi physical properties in the hypothesis of a 
two components model, and we conclude in § [5] 

2. The HI emission 

2.1. GBT observations of the FN1 filament 

The present study focuses on a very diffuse Hi filament lo- 
cated in t he extra-galactic window known as Firback North 
1 (FN1) (jDole et al.l . l200lh located at Galactic coordinates 
(1=84.0°, b=45.1°). The 21 cm Hi emission observations 
used in this analysis were obtained with the Green Bank 
Telescope (GBT). The size of the sub-field of FN1 observed 
is 1.1° x 1.1°, the half-power beam width is 9 '.2 and the 
spectral coverage spans the range from -169 km s^ 1 to 
69 km s^ 1 , with a channel width of 8u = 2.06 km s _1 
and a spectral resolution of 2.5 km s _1 . The noise level in 
a channel is 0.035 K. 

The 21 cm emission of this small patch of the sky is 
characterized by a single spectral component at local veloc- 
ities typical of Hi emission at high Galactic latitude, with a 
narrow core and extended wings (see Fig.[T]). The GBT has 
also detected a faint High- Velocity Cloud (HVC) centered 

at u 120 km s _1 . The average Hi column density of 

the HVC over the field is 0.4 x 10 19 cm~ 2 with a maximum 
column densitjQ of 1.4 x 10 19 cm~ 2 . The 21 cm spectra are 
also characterized by a weak (~ 0.1 — 0.2 K) Intermediate 
Velocity Cloud (IVC) component centered at u ~ —50 km 
s _1 (see Fig. QJbottom) . 

2.2. Properties of the Hi emission spectra 

The 21 cm emission profiles observed are shaped by den- 
sity, velocity and thermal fluctuations of interstellar neutral 
hydrogen on the line of sight. In this section we define quan- 
tities that can be deduced from the observations and that 
will be used further to constrain the physical properties of 
the region observed. 

2.2.1. Brightness temperature of an optically thin line 

In the optically thin limit where Hi self-absorption is neg- 
ligible (|Spitzerl . |1978[ ). the brightness temperature (in K) 
observed at position (x, y) on the sky and in the velocity 
range between u and u + du is: 

1 f H 

T B {x,y,u)du = B{x,y) ® ^ ^ - J n(x,y,z)dz 

1 (u- v(x',y',z)) 2 
X V^ eM 2A* )du - W 

1 The HVC component was computed by integrating over the 
channels between -168.5 and -96.4 km s _1 . 




V (km/s) 

Fig. 1. GBT average 21 cm spectrum of the FN1 field 
shown in linear-linear (top) and linear-log (bottom) scales. 
The spectrum in linear-log scale helps to highlight the faint 

IVC (it 50 km s- 1 ) and HVC (it 120 km s" 1 ) 

components. 



where n(x,y,z) and v(x,y,z) are the three-dimensional 
density and velocity fields. The term A 2 = k T(x, y, z)/m + 
&p S f represents the thermal broadening of the 21 cm line 
and the instrumental broadening; T{x, y, z) is the 3D tem- 
perature field, B(x,y) is the instrumental beam, H is the 
depth of the line of sight, m is the hydrogen atom mass 
and k is the Boltzmann constanfl Even when radiative 
transfer effects are negligible, like for the diffuse filament 
observed here, the 21 cm emission depends in a complex 
manner on the density, velocity and temperature distribu- 
tion on the line of sight. In addition s everal studies showed 
(e.g. jMiville-Deschenes et al.l (|2003aD ) that the Hi density 
and velocity distributions are self-similar, which reflect the 
impact of turbulence on the interstellar medium dynamics. 

One of the challenges of the analysis of 21 cm emission 
is to disentangle density, velocity and temperature fluctua- 
tions, using quantities averaged along the line of sight. To 
constrain the properties of the density, velocity and tem- 
perature of the gas in 3D we use the integrated emission 
map Nu(x,y), the centroid velocity map C(x,y) and the 
velocity dispersion map cr(x, y). 

2.2.2. Integrated emission 

In the optically thin limit, the integrated emission is pro- 
portional to the column density (in cm -2 ), which is simply 
the integral of the density fluctuations on the line of sight: 

N H (x,y) = 1.823 x 10 18 ^ T B (x, y,u)5u (2) 

u 



2 here we consider that the spectral response follows a 
Gaussian function 

3 Here we make the assumption that the depth of the cloud is 
small compared to the distance to the cloud, which allows us to 
do the integral in Cartesian coordinates and not take into ac- 
count the fact that a beam sees structures of increasing physical 
size with distance. 
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Fig. 2. 

FN1. 



Hi column density (left), centroid velocity (middle) and velocity dispersion (right column) of the local Hi gas in 



B(x,y)<g> / n(x' ,y' , z) dz. 
Jo 



(3) 



2.2.3. Velocity centroid 

The velocity centroid map is the average of the velocity as 
weighted by the brightness temperature at each position on 
the sky: 

Y, u uT B (x,y,u) 



C(x,y) 



(4) 



T,u T B(x,y,u) 

In the optically thin limit, the velocity centroid map is equal 
to the following expression: 

„, s B(x,y)® J H n(x' 1 y' 1 z)v{x / ,y' 1 z)dz 

C(x,y) = - . (5) 

B(x,y)®f n(x',y',z)dz 

The velocity centroid map is thus the density-weighted av- 
erage velocity on each line of sight. It depends non- linearly 
on the density and velocity distribution in 3D, but not on 
the gas temperature. 



2.2.4. Velocity dispersion 

The second velocity moment, the velocity dispersion map, 
is the velocity dispersion on the line of sight. It is computed 
the following way: 



<7 2 (x,y) 



J2 U u 2 T B (x,y,u)Su 



C 2 (x,y). 



(6) 



Y JU T B{x,y,u)5u 

This quantity is more affected by instrumental noise but it 
allows one to estimate the gas temperature. It is related to 
the 3D quantities by the following equation: 

<T 2 (x,y) = a 2 urb (x,y) + a* af (x,y) + a 2 herm (x,y) (7) 

where a ps f is the instrumental broadening due to the finite 
spectral resolution, 



B(x, y) g j" n{x', y' , z) v 2 (x', y', z) dz 



2.3. Statistical properties deduced in FN1 

The Hi integrated emission, centroid velocity and veloc- 
ity dispersion maps of the FN1 sub-field, computed respec- 
tively using Eq. [2l Eq. [4] and Eq. [6l are shown in Fig. [2l 
These three maps were computed by using only data points 
with a value higher than 2 times the noise level in a sin- 
gle channel. The average and standard deviations of these 
three mapifl are given in the first row of Table [1] (labeled 
"Total" ) . To estimate the uncertainties listed we have com- 
puted the same three maps by selecting data points with a 
value higher than 1 time and 3 times the noise level. 



2.4. Discussion 

Without any 21 cm absorption measurements in that field, 
our idea was to use the statistical properties of the pro- 
jected quantities together with the spectral shape of the 
21 cm line (Fig. [T]) as the basis for some constraints on 
the level of turbulence. More specifically we wanted to es- 
timate the relative contributions of the temperature and 
turbulence fluctuations to the line profile and to study if 
the observed profiles can be the result of relatively cool gas 
with strong turbulent motions. As stated previously the sta- 
tistical quantities given in the first row of Table [T] represent 
a mixture of density, velocity and temperature fluctuations 
on the line-of-sight. They depend on the statistical proper- 
ties of the 3D density, velocity and temperature fields, and 
of course on the depth of the line-of-sight. 



a turb( x >y) 

B{x,y)®J Q n{x',y',z)dz 
is the broadening from turbulent motions and 



k B(x, y) <g> J Q n{x' , y' , z) T(x' , y' , z) dz 



therm 



(x,y) 



m B{x,y) <£> J Q H n(x' ,y' , z) dz 
is the thermal broadening. 



(9) 



3. Simulation of the GBT observation 

In this section we describe how we produced simulated GBT 
observations to constrain the properties of the 3D density n 
and velocity v fields as well as the average kinetic tempera- 
ture of the gas. One solution to simulate GBT observations 
would be to use hydrodynamics or magneto-hydrodynamics 
simulations of the Hi that would includ e the thermal insta- 
y%&}ty (e.g.. lAudit fc HennebeTIi (|2005fl ). Such simulations 
give a close representation of the actual physics. On the 
other hand they do not provide a direct handle on the power 
spectrum of the density and velocity, nor on the average 
density and velocity, or on the axis ratio of the cloud (i.e., 



4 The standard deviation of a(x, y) is not given here and not 
used in this analysis as it is more affected by noise 
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Table 1. Statistical properties of the local Hi emission. 





<N H > 


<7jV 


< C > 


oc 


< a > 




(10 19 cm" 2 ) 


(10 19 cm" 2 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


Total 


8.1 ±0.1 


0.98 ±0.02 


-5.2 ±0.2 


0.9 ±0.1 


11.3 ±0.6 


Narrow 


2.48 ± 0.01 


0.77 ± 0.03 


-7.27 ±0.01 


0.75 ±0.02 


4.96 ±0.01 


Wide 


5.89 ±0.01 


0.5 ±0.1 


-2.96 ±0.01 


1.66 ±0.06 


13.44 ± 0.01 



The Total row gives the statistical properties computed on the total emission (excluding the HVC). The uncertainties on the Total 
values represents their variance when computed by selecting data points with a 1 and 3 times the noise level threshold. The Narrow 
and Wide rows give the statistical properties of the two components resulting from the Gaussian decomposition (see § I4.1|l . In 
these cases the uncertainties were computed using the uncertainties on the parameters of the fit (see text for details). 



W 21m (Kkm/s) C 2 , cm (km/s) 

15 100 200 300 400-30 -20 -10 




k (deg') k (deg ') 

Fig. 3. Power spectrum analysis of the Hi emission in 
the FN1 region. Top: 21 cm integrated emission (left) 
and centroid velocity map (right) obtained wit h the 
Leiden/Dwingeloo data (jBurton k. Hartmannl . Il99€ ). The 
two maps were computed using velocity channels from - 
86.7 km s _1 to 47.2 km s _1 to avoid contamination by high- 
velocity clouds. Bottom: Power spectrum of the above 
maps with corresponding power law fit. Both power spec- 
tra were divided by the instrumental response of the Leiden 
antenna (beam of 30 arcmin FWHM). 



depth over extent on the sky). Finally doing such simula- 
tions is computationally intensive and it is not practical to 
make a large number of them to estimate statistical vari- 
ance of the results. 

To have explicit control on the statistical properties of 
the fields and the axis ratio of the cloud, and to be able to 
make a large number of realizations, we chose to simulate 
3D velocity and density fields using fractional Brownian 
motion (fBm), which are Gaussian random fields wi th a 
given power spectrum (|Miville-Deschenes et all l2003b[ ) . To 
simulate the GBT observation we made the assumptions 
that we are observing a single Hi cloud on the line of sight, 
that its 3D density and velocity fields are self-similar and 
characterized by power law power spectra and that the gas 



Y | 



L 




Fig. 4. Simulation of a 65 x 65 x 128 density field with 
cr(n)/ < n >=1 and j n — —3.4. The figure shows the 
density fluctuations at surface of three sides of the cube. 



structure and kinematics are statistically isotropic in three 
dimensions. In this framework the statistical properties 
(mean and variance) of the observed quantities (Nh(x,u), 
C{x,y) and a{x,y)) will depend on a relatively limited 
number of parameters: the index of the density and ve- 
locity power spectra (7„ and -f v ), the 3D average density 
< n >, the 3D density dispersion cr n , the 3D average ve- 
locity < v >, the 3D velocity dispersion a v , the average 
gas temperature < T >, the cloud depth along the line of 
sight H and the distance to the cloud D. In fact we reached 
the conclusion that for given value of 7 ra , and to reproduce 
given < Nh > and o~n, there is a close one to one relation 
between the required ratio a n j < n > and the cloud axis 
ratio H/L, where L is the physical size of the field on the 
plane of the sky, that is independent of the distance D to 
the cloud. The model a v and observed a are also indepen- 
dent of distance and only depend on the density contrast 
a n / < n >. It is thus convenient to produce simulations 
for a constant distance D (we adopted 100 pc) but with a 
variable H/L. 
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Fig. 6. Comparison of simulated and observed spectra. Top: Typical simulated spectra, averaged over the 1.1° x 1.1° 
field, for 5 values of the axis ratio H/L. These spectra represent only the turbulent part of the broadening. Bottom: 
Observed average GBT spectrum (diamonds) compared to the average simulated spectra in which thermal and instru- 
mental broadening were added to the turbulence to reproduce the observed < a >. 



3.1. Power spectrum of the Hi emission 

One important ingredient of our analysis is the power 
spectrum of the density and ve locity fluctuations in 3D. 
IMiville-Deschenes et al. (|2003bl ) showed that the power 
spectrum of 21 cm integrated emission and of the cen- 
troid velocity fields can be used to estimate the spec- 
tral in dexes j n and 7„ of the 3D d ensity and velocity 
fieldfl Miville-Desch enes et al.l (|2003a|) used this property 
to estimate the power spectrum of velocity fluctuations 
in 3D of a high-latitude cloud. Such power spectra can- 
not be computed solely with the GBT observations which 
have a very limited spatial dynamic range (from 9.2' to 
1.1°). Instead we used the 2 1 cm Leiden/Dwingeloo data 
(jBurton fc Hartmannl . 11994) of a 40° x 40° region cen- 
tered on the GBT field. The integrated emission and cen- 
troid velocity maps extracted from these data are shown 
in Fig. [3] together with their respective power spectrum. 
We have taken great care here to select velocity channels 
with emission from the local gas only. The power spec- 
tra of the integrated emission and centroid velocity are 
both well fitted by a slope of —3.4 ± 0.2 on scales smaller 
than 5°. At larger scales the power spectra flatten slightly 
which could be partly due to an effect of the tangent 
plane projection. It could also be an indication that we 
have reached the angular scales on the sky that correspond 
to the d epth of the medium (i.e . . the t r ansition from 2D 
to 3 D) (IMiville-Deschenes et all . l2003bl : lElmegreen et all 
[20011 itneaTls et all l2004j T~ 



It is interesting to note at this point that the spec- 
tral indices measured here are compatible with what 
has been measured in anoth er high-latitude field by 
IMiville-Deschenes et al.l (|2003al ) also using 21 cm observa- 
tions. In the GBT simulations we consider that the power 



5 As shown bv lBrunt fe Lowl (|2004f ) and lEscmivel fc Lazarianl 
(2005) this is only true for low contrast density fields (<J„/ < 
nXl). 



spectrum indices of the 3D density and velocity fields are 
In — 7u = — 3.4 following what is observed at larger scales. 

3.2. Simulations of density and velocity fields 

We made simulations of n{x, y, z ) and v{x,y 1 z) using fBm 
(iMiville-Deschenes et all l2003bl ) with j n = j v = -3.4. We 
performed 150 simulations each for 20 different values of 
the H/L ratio (from 0.05 to 10, equally spaced logarithmi- 
cally). The size of our cubes were from 65 x 65 x 3 pixels 
(for H/L = 0.05) to 65 x 65 x 650 pixels (for H/L = 10). 
To minimize the effect of periodicity of fBm objects each 
cube of size 65 x 65 x N was in fact extracted at a ran- 
dom position from a larger fBm cube of 1024 3 pixels. For 
each of these 3000 simulations the n{x, y, z)) and v(x, y, z)) 
fields were constructed independently with no correlation 
between thenfl 

Each pair of (n, v) cubes were then combined to produce 
a brightness temperature cube Tb{x, y, u) (using Eq. [1} on 
the grid of the GBT observation (48 x 48 pixels times 116 
channels of width 6u = 2.06 km s _1 ) and taking into ac- 
count the spatial and spectral instrumental response of the 
GBT. The pixel size in our original 64 x 64 x N cube is 1.1' 
which has to be compared with the GBT beam of 9'. We 
have tested that this was enough spatial resolution to take 
into account fluctuations at scales smaller than the beam. 

To normalize the cubes n and v we computed the col- 
umn density (Nn(x,y)) and centroid velocity (C(x,y)) 
maps from the simulated TB(x,y,u) using equations [5] and 
31 For the velocity, the normalized field is simply 



Av + B. 



(10) 



6 IMiville-Deschenes et al.l j|2003bl ); iBrunt fc Lo^l (|2004f ) 

showed that the level of correlation between the density and 
velocity does not modify significantly the statistics of projected 
quantities. 
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Fig. 5. One component model : physical parameters (den- 
sity contrast a(n)/ < n >, average turbulent broadening 
< oturb >, average temperature < T3 >= T/1000, average 
Mach number < M >) as a function of the axis ratio H/L. 
Also given are the average density < n > and average pres- 
sure < P > for distances to the cloud of 75 pc (solid), 150 
pc (dotted) and 450 pc (dashed). 



We find A and B such that our simulation reproduces the 
observed values < C > and a c . 

For the density field the normalization process is more 
complicated because it has to be positive everywhere. Being 
Gaussian, fBm fields are limited to contrast lower than 
1/3 as 99% of the fluctuations are within 3<r of the mean. 
Interstellar density fields are expected to have higher con- 
trast than 1/3 which calls for a modification of the usual 
fBm. To produce positive and higher-contrast fBm we use 
the following method 



where n m i n is the minimum value of n. We estimate a and 
(3 so that the simulated column density map has the same 
< Nh > and ctjv as the observations. We have tested that 
this modification of the fBm does not impact on its power 
spcctrurrQ. An example of a simulated 3D density field is 
shown in Fig. |4j 

This method is very similar to the exponentiation 
metho d described bv lElmegreenl (|2002T ) and lBrunt &: Heverl 
(2002) but it seems closer to the statistical properties of in- 
terstellar emission. At a given scale, our method produces a 
fluctuations distribution with a larger FWHM and less ex- 
treme values (lower skewness and kurtosis) which is closer 
to the bright ness fluctuations of dust emis sion at 100 /iin 
for instance (jMiville-Deschenes et all l2007t ). 



3.3. Physical properties 

For each normalized simulation we obtain directly a set of 
quantities that are independent of the distance D to the 
cloud: the density contrast (<j n / < n >)@, the standard de- 
viation of the velocity field (cr v ) and the turbulent broad- 
ening on each line of sight (cr tur f ) (a;, y) using Eq. [5]). By 
comparing the average < a 2 urb > to the observed averaged 
line width < a 2 > (both quantities being averaged over the 
field) one can estimate the average gas temperature 



. in/ 
< T >= — (< (7 > 
k 



2 a 2 



psf~ 



< &turb > 



(12) 



and the average Mach number at the scale of the whole 
cloud 



< M >= a v 



k<T> 



(13) 



where a v is the standard deviation of the 3D velocity field 
v . These are also independent of D. 

Assuming a distance D to the cloud, L = Dtwa.9 where 
9 is the angular size of the GBT field (1-1°) one can also 
compute the average gas density 



< n >= 



L < N H > 



(14) 



H Dtim6 

and the average pressure < P >=< n >< T >. 
3.4. Results 

The results of these normalized simulations are compiled 
in Fig. On the top four panels are the density contrast 
a n j < n >, the turbulent broadening aturbi the average 
temperature < T > - in units of 1000 K - and the aver- 
age Mach number < M >, as a function of the axis ratio 
H/L. The solid line is from the average for the 150 realiza- 
tions done at a given H/L and the error bars represent the 
standard deviation. 

The bottom two panels of Fig. Ogive the average density 
and pressure as a function of H/L and for three distances 
D = 75, 150 and 450 pc. In the absence of any measurement 
of the distance to the cloud we have estimated a range and 



^' = a{n — n 



HI i 11 I 



(11) 



7 From our tests on cubes with 0.1 < H/L < 10 and a spectral 
index 7 = —3.4 we did not observe significant changes of the 
power spectrum after the modification of the fBm using Eq. 1111 

8 Both the average density < n > and the standard deviation 
of the density cr„ depend on D but not their ratio 
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Fig. 8. Column density (left column), centroid velocity (middle column) and velocity dispersion (right column) of the 
21 cm emission of the narrow (bottom) and wide (top) components in the local-velocity gas. Note the widely different 
ranges indicated on the color bars. 



a most probable value. lLockman fc Gehmanl (|l99ll ) showed 
that the scale height (HWHM) of Hi is 125 pc in the solar 
neighborhood. At the Galactic latitude of FN1 (b = 45.1°) 
we can consider that most of the clouds lies within 150 pc 
(la of the scale height) from the Sun, and that 99% of the 
Hi is closer than 450 pc. In add ition, using a ma pping of 
Nal absorption against local stars lSfeir et al.1 (|1999D showed 
that there is very little neutral gas at a distance < 75 pc 
from the Sun. So we consider 75 and 450 pc as the lower and 
upper limit on distance and 150 pc as the most probable 
value. 

As expected the required density contrast and the tur- 
bulent broadening increases with H/L. The brightness fluc- 
tuations average out as the depth of the line of sight in- 
creases, and so stronger density contrasts are needed to 
explain the observed level of column density fluctuations 
((7jv)- The same is true for the turbulent broadening: for 
deeper clouds the central limit theorem tends to wash out 
the differences of centroid velocity between two adjacent 
lines of sight, and therefore, to explain a given level of 
fluctuations in the centroid velocity map (<tc)j one needs 
stronger turbulent motions. 

3.5. Line profiles 

Another important way of analyzing the results of these 
simulations is to inspect the shape of simulated 21 cm spec- 
tra. In Fig. [5] we show typical simulated 21 cm spectra, 
averaged over the 1.1° x 1.1° field, for different values of 
H/L. On the top row the spectra represent what would be 
observed with no thermal and instrumental broadenings: 
the spectral structure is only due to turbulent motions. As 
discussed, with the increase of the depth of the cloud the 
spectra get wider but also smoother due to the averaging 



of density and velocity fluctuations on the line of sight. 
Compared to the GBT observations in the bottom row, this 
turbulent line profile is very narrow. For this specific field, 
we have to reject the hypothesis that the observed spectra 
is a product of highly turbulent cold gas: the fluctuations 
seen in the centroid velocity maps are to small for that. 

On the bottom row, we show the simulated spectra 
when thermal and instrumental broadening are added (solid 
line) to reproduce the observed line width; these curves 
have the same velocity dispersion, as defined by Eq. [6l as 
the average GBT spectrum (diamond). The most important 
result here is the fact that the simulated spectra do not re- 
produce well the observations. The simulated spectra in the 
bottom row are very close to Gaussian: the turbulent con- 
tribution being relatively small, strong thermal broadening 
(T = 13000 - 15000 K) is needed to get to the a observed. 
More specifically the simulated spectra failed to reproduce 
the peaked central part and the extended wings of the ob- 
served spectrum. The conclusion drawn from this analysis 
is that gas at a single hot temperature cannot reproduce the 
observed spectra adequately either. The fact that the ob- 
served spectra have a relatively narrow core and extended 
wings is instead suggestive of two components, one colder 
than the other. In the following we explore this possibility. 

4. Two components model 

4.1. Gaussian decomposition 

Many studies of 21 cm observations rely on a Gaussian 
decomposition of the spectra and make the assumption 
that each Gaussian componen t is a physical structure (a 
cloud) on the line of sight fe.g.. lHeiles fc Trolandl (|2003al )V 
On the other hand, it is not always conceptually easy to 
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Fig. 9. Uncertainty maps of the maps shown in Fig. [51 




-60 -30 30 60 -60 -30 30 60 -60 -30 30 60 

v (km/s) v (km/s) v (km/s) 



Fig. 7. Example of the Gaussian decomposition of the 21 
cm spectra. The Hi integrated emission is shown in the cen- 
tral panel, surrounded by 8 typical spectra (black line) with 
their associated Gaussian decomposition (red lines). Three 
components were used to fit the data : two local components 
(one narrow and one wide) and one intermediate velocity 
component. 



reconcile such a decomposition with the fact that the in- 
terstellar medium has structure at all scales. In fact the 
Gaussian decomposition is facilitated by the significant 



thermal broadening of the 21 cm line and by the instru- 
mental response (spectral and spatial) which often have 
quasi-Gaussian functions. In addition, with an increase of 
the path length along a line of sight, the central limit theo- 
rem also tends to bring the velocity and density fluctuations 

closer to a Gau ssian distribution. 

As stated bv lHeiles fc Trolandl (|2003bh the Gaussian fit- 
ting procedure is also subjective and non-unique. To make a 
Gaussian decomposition of the data we used a constrained 
fit method that 1) limits the parameter space (centroid and 
sigma) for each component and 2) iterates to find a spa- 
tially smooth solution for N, < v > and a, which reduces 
slightly the subjectivity of the fit and certainly increases its 
robustness. 

The spectra are in fact very well fitted by a sum of 
three Gaussian components, a wide one typical of the 
Warm Neutral Medium ( WNM) , a narrow one representing 
colder gas and a very weak one at intermediate velocity. 
Typical fits are preseted in Fig. [7J The intermediate ve- 
locity component is very weak (average column density is 
0.8 x 10 19 cm~ 2 , maximum is 1.3 x 10 19 cm -2 ). It probably 
does represent a physically distinct component and will not 
be used in the following analysis. The resulting column den- 
sity, velocity centroid and velocity dispersion found for the 
two local components are shown in Fig. [8j The uncertainty 
maps computed using the uncertainties on the fit parame- 
ters are shown for each quantity in Fig. [9l The statistical 
properties of the integrated emission, centroid velocity and 
velocity dispersion of the narrow and wide components are 
given in Table [II The uncertainty on these statistical prop- 
erties were computed using 1000 realizations of each param- 
eter maps shown in Fig.[8jto which we have added a random 
value at each position. The random value was taken from a 
Gaussian distribution of width given by the uncertainty on 
the fit parameters given in Fig. [9] The uncertainties given 
in Table [T]are the variance (for the average) or the bias (for 
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Fig. 10. Two components model : physical parameters (density contrast cr(n)/ < n >, average turbulent broadening 
< Cturb >j average temperature < T3 >= T/1000, average Mach number < M >) as a function of the axis ratio H/L. 
Are also given the average density < n > and average pressure < P > for distances to the cloud of 75 pc (solid), 150 pc 
(dotted) and 450 pc (dashed). 



the standard deviation) of the resulting statistical proper- 
ties. 

This decomposition shows that the column density of 
the field is dominated by the wide component (~ 2/3). On 
the other hand the fluctuations of the total integrated emis- 
sion are dominated by the narrow component. For the dis- 



cussion in § [3] this is compatible with the wide component 
having a depth much larger than the narrow component. 
We note a significant increase of the wide component col- 
umn density at the tip of the main filament of the field, 
which does not show up in the narrow component column 
density. On the other hand the narrow component reaches 
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its smallest velocity dispersion (2.5 km s" 1 ) at that position 
which coincide spatially with an interesting convergence of 
the centroid velocity of the warm component. This spatial 
correlation of the column density of the wide component 
and the velocity dispersion of the narrow component at the 
tip of the filament is significant and can't be explained by 
an artifact of the fitting procedure. As it can be seen in 
Fig. [HI this is the location where the constraints on the fit 
are the strongest. 

4.2. Simulations 

To estimate the physical properties of the narrow and wide 
components we used the exact same method as described 
previously. We made fBm simulations of 3D density and 
velocity fields in order to reproduce the observed statis- 
tical properties of the narrow and wide components, as a 
function the axis ratio H/L. In the absence of constraints 
we made the assumption that the narrow and wide com- 
ponents are both characterized by the same power spec- 
trum (7„ = 7„ = —3.4) and that they are uncorrelated. 
Like for the previous simulations we made 150 realizations 
per H/L value both for the narrow and wide components. 
We also made 30 realizations for j n = j v = —3.0 and for 
7n = 7« = —3.8 in order to test the robustness of our results 
on the uncertainty on the power spectrum slope. 

4.3. Results 

The results of the simulations are compiled in figures ITD1 
and llll Like for the simulations done in the one-component 
hypothesis, the density contrast and turbulent broadening 
increase with H/L but with different average values. At 
a given H/L the narrow component has a higher contrast 
and a lower turbulent velocity dispersion than the wide 
component. 

The effect of the uncertainty on the power spectrum 
slope on these two quantities is shown in 1111 There is a 
bias introduced by the value of 7„ and j v used. At a given 
H/L, steeper (flatter) power spectra leads to lower (larger) 
density contrast and turbulent broadening. Nevertheless, 
even for the extreme values of the power spectrum index 
used here (-3.0 and -3.8) we find very similar results than 
for 7„ = j v — -3.4. 

Typical simulated spectra are compared to the average 
observed spectrum in Fig. [T2] As in Fig. O on that latter 
figure the top panels show typical narrow (blue) and wide 
(red) spectra if there was only turbulence contributing to 
the broadening. On the bottom panels thermal and instru- 
mental broadening were added to the narrow and wide spec- 
tra to match the corresponding velocity dispersions < a > 
given in Table [TJ Contrary to the one component model 
(see Fig.[S]) the sum of the two components, shown in black, 
matches very well the observations, shown as diamonds. 

This simulation, based on a simple representation of 
the density and velocity fields (self-similar), gives insight 
into the physical parameters of the Hi in this very dif- 
fuse region of the Galaxy. If we consider a typical pressure 
(P ~ 5000 K cm" 3 ) and distance (D ~ 150 pc) for this 
high latitude cloud, the bottom panels of Fig. [TU] suggests 
that the H/L ratio of the narrow component is about 1 but 
reaches 10 for the wide component. The wide component 
is therefore much more spread out on the line of sight than 



1.0 




™i — ■ ■ i — i_r " 

One component 






0.1 



1.0 10.0 0.1 
H/L 



1.0 

H/L 



10.0 0.1 




10.0 



Fig. 11. Effect of a variation of the power spectrum in- 
dex (7„ and 7„) for the one component model (left col- 
umn) and the two components models (middle and right 
columns). Result of the simulation for the density contrast 
<j(n)l < v > and turbulent broadening < a tur b > when 
one considers a power spectrum index of j n — -f v — —3.0 
(dashed) and 7 n = 7-u = —3.8 (dotted) (30 realizations per 
H/L value). The solid line with error bars is the results for 
7„ = 7„ = —3.4 (150 realizations per H/L value - the same 
curves as shown in figures [5] and I10p . 



the narrow one which is more confined spatially. In ^ 13.11 wc 
noted a break in the Leiden/Dwingeloo data power spec- 
trum at a scale ~ 10° (see Fig. [3]), also an indication that 
H/L ~ 10 for the wide component which dominates the 
column density. 

Using these values of H/L as basis, we found that 
a n / < n >~ 0.8 for both phases, which are relatively low 
contrast valuefl These values of density contrasts are com- 
patible with what is u sually obtained in M HD simulations 
of interstellar clouds (|Brunt fc Lowi . l2004h . We also find 
that the turbulent broadening is ~ 2 km s _1 for the nar- 
row component and ~ 10 km s _1 for the wide component. 

In the hypothesis that H/L = 1 for the narrow com- 
ponent and H/L = 10 for the wide component, the wide 
component has a temperature (T ~ 10000 K) and a density 
(< n >~ 0.5 — 1) which are typical of WNM gas. On the 
other hand the narrow component has an average temper- 
ature (T ~ 2300 — 2800 K) that is clearly not representa- 
tiv e of CNM gas (th e media n CNM temperature obtained 
bv lHeiles fc Trolandl (|2003bh is 70 K). The temperature of 
the narrow component is in fact in the thermally unstable 
range. Globally CNM type gas does not contribute a large 
fraction of the column d ensity in the GBT field. This is 
not totally unexpected as iHeiles Ez Trolandl (|2003b| ) found 
several lines of sight with no CNM. This analysis in fact 
leads to the conclusion that, in this field, ~ 30% of the to- 
tal Hi column density is in the thermall y unstable regime. 
This is similar to the average obtained bv lHeiles &; Trolandl 
(|2003b| ) using a completely different technique. In addition, 
we argue that the unstable phase has a lower filling factor 



These relatively low contrasts give an indication that the 
velocity power spectrum obtained from the centroid velocity is 
not affected significantly by density fluctuations. 
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Fig. 12. Comparison of simulated and observed spectra. Same as Fig. [6] but for the narrow (blue curve) and wide (red 
curve) components. The black curves on all plots is the sum of the narrow and wide components. 



(lower H/L) than the WNM. In that respect the unstable 
gas would have a structure closer to the CNM. 

Regarding the turbulent motions, the unstable gas is 
definitely sub-sonic with < M >^ 0.3. On the other hand 
the WNM gas is quite turbulent and likely to be supersonic 
with < M >~ 1. This could give an indication on the 
dissipation of turbulent motions associated with the phase 
transition (WNM^unstable^CNM). 

Globally our results are in good ac cord with the simu- 
lations by lAudit & Hennebelld (|2005l ) of the Hi thermal 
instability. They show that, for slightly turbulent fields, 
up to 30% of the gas is in the thermally unstable phase 
(200 < T < 5000 K) with an average density n — 5 cm" 3 . 
For H/L — 1 and D ~ 150 pc we also find < n >~ 
5 c m" 3 for our narrow compon ent. It is also interesting 
that lAudit k Hennebelld (|2005f) found the thermally un- 
stable gas to be filamentary in accord with the low H/L 
ratio of the observed narrow component. 



4.4. Formation of a CNM structure ? 

The analysis above gives average values over the field of 
view but one can also estimate the physical parameters of 
a specific structure seen in the column density map, like 
the main filament for instance which has a low velocity dis- 
persion and could contain classical CNM gas. This filament 
has an angular transverse size of 10 arcmin (which is barely 
resolved by the GBT). 

Considering the properties of the narrow component at 
the position of minimum a (Nh ~ 3.3 x 10 19 cm -2 and 
a = 2.5 km s" 1 ) and assuming H/L = 1 one can derive its 
depth 



^ = 0.44^50 (pc), 
its average density 

< n >= 24/£>i 50 (cm -3 ) 
average temperature 

<T>= 204^000 ^150 (K) 



(15) 
(16) 
(17) 



turbulent broadening 

aturb = a/5-1 - 1.7P 50 oo D 150 (km/s) 
and average Mach number 



< M >= 



-P5000 L>i 



- 1 



no 



(18) 



(19) 



where D 150 = D/(150 pc) and P5000 = P/(5000 K cm" 3 ). 
For those typical values of distance and pressure we get 
Pturb = 1-8 km s" 1 and < M >= 1.4. These physical pa- 
rameters are typical of CNM gas, we therefore conclude that 
it is likely that a significant fraction of the column density 
at the tip of the filament seen in the FN1 field is CNM 
gas. The fact that we observe opposite velocity gradients of 
WNM gas at the same position (see Fig. [51) is reminiscent 
of wh at is seen in the simulations of lAudit fc Hennebelld 
(2005|) where CNM structures are efficiently formed in a 
large scale convergent flow. 



5. Conclusion 

In this paper we presented a new technique to constrain the 
physical properties of the 3D density and velocity fields of 
Hi clouds. We showed that from a set of statistical quanti- 
ties computed with the observations (< Nh >, o^r, < C >, 
crc, < cr >, 7jv, 7c) and the spectral shape of the 21 cm 
line, it is possible to constrain the physical properties of 
the Hi on the line of sight like density contrast, turbulent 
broadening, average temperature and average Mach num- 
ber. 

The basis of the method is to constrain the underly- 
ing physical properties of the cloud by producing a large 
number of 3D realizations of the density and velocity fields 
with proper statistical properties. This type of method is 
routinely used in the analysis of the cosmic microwave back- 
ground (CMB). One important difference with the CMB 
analysis is the fact that we do not have a complete descrip- 
tion of the statistical properties of interstellar fields and 



12 



Miville-Deschenes, M.-A. & Martin, P.G.: A diffuse HI structure 



that Gaussian random fields (or fBm) are not a good repre- 
sentation of the observations. One known limitation of fBms 
for the analysis of interstellar matter is the fact that they 
are isotropic and Gaussian objects. In this analysis we used 
a method to modify classical fBm objects such that they re- 
produce the power spectrum and the non-Gaussian proper- 
ties of density fields observed in numerical simulations and 
observations of interstellar clouds. A further step would be 
to produce such artificial fields with anisotropic structure 
(obvious in interstellar emission) and a better description 
of the non-Gaussianity of the velocity field (weaker than 
the non-Gaussianity of the density field). Nevertheless we 
argue that the statistical properties used here to estimate 
the 3D physical properties of the gas are more controlled by 
the power-law power spectrum of the density and velocity 
fields (and by the non-Gaussian properties of the density 
field) rather than by their anisotropy. 

We applied this new technique to analyze 21 cm obser- 
vations obtained at the GBT of the FN1 field, a very dif- 
fuse region at high Galactic latitude used for extra-galactic 
studies. First we showed that the observed spectra of FN1 
cannot be explained by CNM type gas affected by strong 
turbulent motions. We were able to exclude this hypothe- 
sis by showing that it is incompatible with the statistical 
properties of the centroid velocity. More importantly we 
showed that the shape of the spectra (narrow core and ex- 
tended wings) cannot be well reproduced with only one Hi 
component, whatever its temperature. 

Instead the observations are well fitted by a combination 
of two spectral components: one narrow (< a >= 5.0 km 
s _1 ) that traces thermally unstable gas, and one wide 
(< a >= 13.4 km s -1 ) which has the physical proper- 
ties of the WNM. The results of our analysis show that 
there is no significant CNM type gas in the field and that 
the thermally unstable phase contributes ~ 30% of the col- 
umn density, ne ar the average reported in recent obser- 
vational res ults (|Heiles fc Trolandl l2003bl) and numerical 
simulations (| Audit fc Hennebelld . l2005f) . Both components 
have a density contrast a n / < n > close to 0.8 but the 
WNM component is much more spread out on the line of 
sight with a filling factor 10 times higher than the unstable 
gas. We also conclude that the WNM gas is mildly super- 
sonic (< M >^ 1) and the unstable phase is definitely 
sub-sonic (< M >~ 0.3). 

Finally this portion of the FN1 field contains a narrow 
filament (axis ratio > 4), almost unresolved by the GBT, 
which has a low velocity dispersion (2.5 km s" 1 ). Our anal- 
ysis indicates that this structure, which could contain CNM 
type gas, is found at the locus of opposite velocity gradients 
in the WNM which could have triggered a phase transition. 
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